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Abstract 

We demonstrate that stochastic differential equations (SDEs) driven by fractional 
Brownian motion with Hurst parameter H > ^ have similar ergodic properties as 
SDEs driven by standard Brownian motion. The focus in this article is on hypoelliptic 
systems satisfying Hormander's condition. We show that such systems enjoy a suitable 
version of the strong Feller property and we conclude that under a standard control- 
lability condition they admit a unique stationary solution that is physical in the sense 
that it does not "look into the future". 

The main technical result required for the analysis is a bound on the moments of the 
inverse of the Malliavin covariance matrix, conditional on the past of the driving noise. 

Resume 

Nous demontrons que les equations differentielles stochastiques (EDS) conduites par 
des mouvements Browniens fractionnaires a parametre de Hurst H > ^ ont des pro- 
prietes ergodiques similaires aux EDS usuelles conduites par des mouvements Brown- 
iens. L'interet principal du present article est de pouvoir traiter egalement des systemes 
hypoelliptiques satisfaisant la condition de Hormander. Nous montrons qu'une ver- 
sion adequate de la propriete de Feller forte est verifiee par de tels systemes et nous en 
deduisons que, sous une propriete de controUabilite usuelle, ils admettent une unique 
solution stationnaire qui ait un sens physique. 

U ingredient principal de notre analyse est une borne superieure sur les moments 
inverses de la matrice de Malliavin associee, conditionnee au passe du bruit. 
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Introduction 

this paper we study the ergodic properties of stochastic differential equations (SDEs) 
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driven by independent fractional Brownian motions (fBm) with fixed Hurst index 
H € (1/2, 1). Recall that fBm is the centred Gaussian process with Bq = and 

E\Bl-B'f = \t-s\^" . 

Although the basic theory of SDEs driven by fBm (for H e (1/4, 1)) is now well es- 
tablished (see for example [Nua06, Mis08, CQ02, FrilO]), the ergodicity of such SDEs 
has only been studied recently. The main difficulty is that non-overlapping increments 
of fBm are not independent, so that we are dealing with processes that do not have 
the Markov property. As a consequence, the traditional ergodic theory for Markov 
processes as in [MT93] does not apply to this situation. 

To our knowledge, the first result on the ergodicity of (1.1) was given in [Hai05], 
where the author considers the case where the vector fields Vi for i > are constant 
and non-degenerate in the sense that the corresponding matrix is of rank n. In this case, 
an explicit coupling argument allows to show that any solution converges to the unique 
stationary solution at speed bounded above by approximately t^^^^. The argument 
however relies in an essential way on the additivity of the noise. 

A more general construction was given in [HO07, Hai09], where the authors built a 
general framework of stochastic dynamical systems (SDSs) very similar to that of ran- 
dom dynamical systems (RDSs) [Arn98]. In this framework, the notion of an "invariant 
measure" can be defined for an autonomous equation driven by an ergodic process with 
stationary but not necessarily independent increments. Besides technical questions of 
continuity, the main difference between the framework of SDSs and that of RDSs is one 
of viewpoint. An RDS is considered as a dynamical system on the state space times the 
space of all 'futures' of the driving process, endowed with a cocycle structure. An SDS 
on the other hand is considered as a Markov process on the state space times the space 
of all 'pasts' of the driving process, endowed with a cocycle structure. This shift of 
viewpoint has two important advantages when studying the ergodic properties of such 
systems: 

• Since the 'future' of the driving noise is not part of the augmented phase space, 
we do not have to deal with 'unphysical' stationary solutions, where the value 
of the process at a given time is a function of the future realisation of the driving 
noise. 

• It allows to build our intuition based on the theory of Markov processes, rather 
than the theory of dynamical systems. There, many criteria for the uniqueness 
of an invariant measure are known to exist. 

Inspired by the celebrated Doob-Khasminskii theorem (see for example [DPZ96]), 
[HO07, Hai09] introduced a notion of a strong Feller property for SDSs that is a nat- 
ural generalisation of the same notion for Markov semigroups. It turns out that for a 
large class of 'quasi-Markovian' SDSs (this includes SDEs of the type (1.1), but also 
many other examples), uniqueness of the invariant measure then follows like in the 
Markovian case from the strong Feller property, combined with some type of topolog- 
ical irreducibility. To illustrate these results, denote by At.x the closure of the set of all 
points that are accessible at time t for solutions to (1.1) starting at x for some smooth 
control B. Then, one has 

Theorem 1.1 If there exists t > Q such that (1.1) is strong Feller at time t and such 
that PlajeM" ■^t,x 7^ 4>, then (1.1) can have at most one invariant measure. 
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Remark 1.2 Since (1.1) does not determine a Markov process, it is not clear a priori 
what is the correct notion of an invariant measure. The notion retained here is the one 
introduced in [HO07, Hai09] . Essentially, an invariant measure for ( 1 . 1 ) is the law of a 
stationary solution which does not 'look into the future' . Similarly, it is not clear what 
is the correct notion of a 'strong Feller property' for such a system. Precise definitions 
can be found in Section 2.1 below. 

Furthermore, it was shown in [HO07] that if (1.1) is 'elliptic' in the sense that 
the matrix determined by {Vi{x)}i^o is of full rank at every x E R", then it does 
indeed satisfy the strong Feller property. In the case of SDEs driven by white noise, it 
is known on the other hand that the strong Feller property holds whenever the vector 
fields Vi satisfy Hormander's bracket condition, namely that the Lie algebra generated 
by {dt + Vq, Vi, . . . , Vd} spans R"+^ (the additional coordinate corresponds to the dt 
direction) at every point. See also Assumption 4. 1 below for a precise formulation. 

Note at this stage that the strong Feller property of (1.1) is closely related to the 
existence of densities for the laws of its solutions, but cannot be deduced from it di- 
rectly. The existence of smooth densities under an ellipticity assumption was shown in 
[HN07, NS06], while it was shown under the assumption of Hormander's bracket con- 
dition in [BH07]. In this paper, we are going to address the question of showing that 
the strong Feller property (in the sense given in [HO07, Hai09], see also Definition 2. 1 
below) holds for (1.1) under Hormander's bracket condition. 

Our main result is: 

Theorem 1.3 Assutne that the vector fields {Vi} have bounded derivatives of all orders 
for i > and furthermore are bounded for i > 1. Then, if Hormander's bracket 
condition holds at every point, (1.1) is strong Feller In particular, if there exists t > 
such that n2:effi" ■'^t.x 7^ ond if there exist constants Mi, M2 > such that 

{x,Vo(x)) < Ah ^ M2\x\\ (1.2) 

then (1.1) admits exactly one invariant measure. 

Remark 1.4 Since we would like to satisfy (1.2) in order to guarantee the existence of 
an invariant measure by [HO07, Prop 4.6], we do not impose that Vq is bounded. This 
creates some technical difficulties in obtaining a priori bounds on the solutions which 
are not present in [BH07]. 

Remark 1.5 The notion of an 'invariant measure' studied in [BC07] is much more 
restrictive than the one studied here, since they considered the situation where the initial 
condition of the equations is independent of the driving noise, whereas we consider the 
case where interdependence with the past of the driving noise is allowed. In particular, 
even the case n = d = 1 with Vo{x) = —x and Vi(x) = 1 admits no invariant measure 
in the sense of [BC07]. 

On the other hand, our definition of an invariant measure is more restrictive than 
that of a random invariant measure in [Arn98]. Since it generalises the notion of an 
invariant measure for a Markov process, this is to be expected. Indeed, there are sim- 
ple examples of elliptic diffusions on the circle whose Markov semigroup admits a 
unique invariant measure, but that admit more than one random invariant measure when 
viewed as RDS's. 
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The remainder of the article is structured as follows. In Section 2, we set up our 
notations and recall the relevant results from [Hai05] and [HO07]. In Section 3, we 
derive the necessary moment estimates for the solutions of the SDE (1.1). In Section 4, 
we then obtain an invertibility result on the conditioned Malliavin covariance matrix of 
the process (1.1), provided that Hormander's condition holds. Similarly to [BH07], this 
allows us to show the smoothness of the laws of solutions to (1.1), conditional on the 
past of the driving noise. In Section 5, we finally show that the strong Feller property 
is satisfied for (1.1) under Hormander's condition, thereby concluding the proof of 
Theorem 1.3. 
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2 Preliminaries 

In this section, we describe the general framework in which we view solutions to (1.1) 
and we recall some of the basic results from [Hai05, HO07, Hai09]. 

2.1 General framework 

Let C^(R_;M'^) be the set of smooth compactly supported functions on ]R_ which 
take the value at the origin. For 7, 5 S (0, 1) define Wi^.s) to be the completion of 
Co°(M_ ; M'') with respect to the norm 

II II Ht)-ui,s)\ 

sdm^ K-sr(i + \t\ + \s\y 

We write yV(7,5) for the corresponding space containing functions on R+ instead of 
K_. Note that when restricted to a compact time interval, the norm (2.1) is equivalent 
to the usual Holder norm with exponent 7. Moreover, is a separable Banach 

space. 

¥oxH e (1/2, 1), 7 e (1/2, iJ)and7+(5 G {H, 1), it can be shown that there exists 
a Borel probability measure P// on = W{^^s) x VV(7,<5) such that the canonical process 
associated to Yh is a two-sided fractional Brownian motion with Hurst parameter H 
[Hai05, HO07]. We fix such values 7, 5, and H once and for all and we drop the 
subscripts (7, 5) from the spaces W and YV from now on. Note that is a product 
measure if and only if iJ = 1/2, but in general we can disintegrate it into a projected 
measure (denoted again by P/j) on W and regular conditional probabilities 7^(w, •) 
on W. Since P/j is Gaussian, there is a unique choice of V{lo, •) which is weakly 
continuous in uj. 

We denote by <^ : R" x >V -)■ C(]R+, M") the solution map to (1.1). It is well-known 
[ You36, HN07] that if the vector fields Vi have bounded derivatives, then this solution is 
well-defined in a pathwise sense (integrals are simply Riemann-Stieltjes integrals) for 
all times. We also define the shift map 6*4: W y.W ^ W y.W hy identifying elements 
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(u! , ijj) £ yV X yV with the corresponding 'concatenated' function U (I>) : M M"^ 
and setting 

etiuj, w) - (cj U u})(t + •) - U w)(-) . 

Denote (ptix,Cj) (p{x,ui){t) and denote by $t : M" x W x VV R" x W the 
augmented solution map given by 

^t(x, uj, ui) = {iptix, w), UwOtiuj, oj)) , 

where Hyv is the projection onto the W-component. For a measurable function / : 
X 1-^ y between two measurable space X, y and a measure fi on X we define the 
push forward measure ~ fi o f~^. With this notation, we can view the solutions 
to (1.1) as a Markov process on M" x W with transition probabilities given by 

Qt{x,uj;-) = ^t{x,uj,-rV{iJ,-). 

These transition probabilities can actually be shown to be Feller [HO07], but they are 
certainly not strong Feller in the usual sense, since transition probabilities starting from 
different instances of ut remain mutually singular for all times. Instead, we will use a 
notion of strong Feller property that is better suited for the particular structure of the 
problem at hand, see Definition 2.1 below. 

However, the question of uniqueness of the invariant measure for (1.1) should not 
be interpreted as the question of uniqueness of the invariant measure for Qt. This is 
because one might imagine that the augmented phase space M" x W contains some 
'redundant' randomness that is not used to describe the stationary solutions to (1.1). 
(This would be the case for example if the V^'s are not always linearly independent.) 
One would like therefore to have a concept of uniqueness for the invariant measure that 
is independent of the particular description of the driving noise. 

To this end, we introduce the Markov transition kernel Q from M" x >V to C(M+ , R") 
given by 

Q(x,uj;-) = ip(x,-rV(uj,-)- 

This is the conditional law of the solution to (1.1) given a realisation ui of the past of 
the driving noise. The action of the Markov transition kernel Q on a measure fj, on 
R" X W is given by 

Qli (A) = / Qix, uj; A) fi{dx, duj) . 

With this notation, we have a natural equivalence relation between measures on R" x W 
given by 

With these definitions at hand, the statement 'the invariant measure for (1.1) is unique' 
should be interpreted as 'the Markov semigroup Qt has a unique invariant measure, 
modulo the equivalence relation ~'. 

2.2 Ergodicity of SDEs driven by fBm 

We now summarise some of the relevant results from [Hai05, HO07, Hai09] giving 
conditions for the uniqueness of the invariant measure for (1.1). This requires a notion 
of 'strong Feller' property for (1.1). We stress again that the definition given here has 
nothing to do with the strong Feller property of Qt- It rather generalises the notion of 
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the strong Feller property for the Markov process associated to (1.1) in the case where 
the driving noise is white in time. 

Let Rt : C(M+, X) i-^ C{[t, oo), X) denote the natural restriction map and let || • Hxy 
denote the total variation norm. We then say that: 

Definition 2.1 The solutions to ( 1 . 1 ) are said to be strong Feller at time t if there exists 
a jointly continuous function £ : (R")^ x W ^ M_|_ such that 

\\R*Q(x,u:;-)- R;Q(y,Lo;-)hv < e(x,y,uj) , (2.2) 

and £{x, x, w) = for every x G M" and every uj G W. 

We also introduce the following notion of irreducibility: 

Definition 2.2 The solutions to (1.1) are said to be topologically irreducible at time 
t if, for every x e M", a; e W and every non-empty open set U C K", one has 
Qt(a;,w;[/ X W) > 0. 

The following result from [HO07, Thm 3.10] is the main abstract uniqueness result 
on which this article is based: 

Theorem 2.3 If there exist times s > and t > Q such that the solutions to (1.1) 
are strong Feller at time t and irreducible at time s, then (1.1) can have at most one 
invariant measure. 

Therefore, in order to prove Theorem 1.3, the only missing ingredient that we need to 
estabUsh is the strong Feller property for (1.1) under Hormander's bracket condition. 

2.3 Notations and definitions 

For T > and measurable / : [0, T] ^ M", set 

||/||o.T.oo= sup \f{t)\, ||/||o,T,^= sup Mi:^. (2.3) 

For a £ (0, 1), we also define the fractional integration operator I" and the corre- 
sponding fractional differentiation operator 2?" by 

Z"/(i)= -i- / (t-s)"-l/(5)ds. 



r(«) Jo 



t (2.4) 

r(l - a) dt Jq 



Remark 2.4 The operators T" and 2?" are inverses of each other, see [SKM93] for a 
survey of fractional integral operators. 

The reason why these operators are crucial for our analysis is that the Markov 
transition kernel V appearing in Section 2.1 is given by translates of the image of 
Wiener measure under I^^^/^, see also Lemma 2.7 below. 

Indeed, by the celebrated Mandelbrot- Van Ness representation of the fBm [MVN68], 
we may express the two-sided fBm B with Hurst parameter H G (0, 1) in terms of a 
two-sided standard Brownian motion W as 

Bt = aH / i-r)"-^/HdWr+t - dWr) (2.5) 

— OO 
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for some an > 0. The advantage of this representation is that it is invariant under 
time-shifts, so that it is natural for the study of ergodic properties, see [ST94] for more 
details. ^ 
Define now the operator Q: W ^ W hy: 

goj{t) = ^H / -g(- W(-r)dr, (2.6) 
Jo r \rJ 

where the kernel g is given by 

g{x) ^ x^-'l^ +(H- 3/2) X f + , (2.7) 

and the constant is given by = {H — l/2)aHCti-H- Here, an is the constant 
appearing in (2.5). It was shown in [HO07, Prop A. 2] that Q is indeed a bounded linear 
operator from W into W. Furthermore, we can quantify its behaviour near t = in the 
following way: 

Lemma 2.5 On any time interval bounded away from 0, the map t h- !■ Gio(t) is C°°. 
Furthermore, if we set fuiit) ~ t-^Quj{t), then we have fuiO) ~ and for every T > 0, 
there exists a constant Mt such that \\fu)\\o,T,-i < -^^tH'^I 1(7,(5). 

The proof of this result is postponed to the appendix. 

Remark 2.6 Since is Gaussian, it follows in particular that both ||tJcj||o,T,7 and 
1 1 1 1 o,T,7 have exponential moments under f'n by Fernique's theorem [Bog98] . 

Let T/, : w ^ w + h denote the translation map on W. We cite the following result 
from [HO07, Lemma 4.2]: 

Lemma 2.7 The regular conditional probabilities V(uj, ■) of¥ff given G W are 

where W is the standard d-dimensional Wiener measure over R_|- and I" is as defined 
in (2.4). 

We will henceforth interpret the above lemma in the following way. The driving 
fractional Brownian motion B can be written as the sum of two independent processes: 

nt 

Bt = Bt + (guj)(t) = Bt+mt, Bt = aH (t- s)^-!/^ dW(s) , (2.8) 

Jo 

where W is a standard Wiener process independent of the 'past' cj e W. This notation 
will be used repeatedly in the sequel. 

Remark 2.8 In the notation of Section 2. 1 , we will repeatedly use the notations E (and 
P) for conditional expectations (and probabilities) over a) with ui fixed. In the notation 
(2.8), this is the same as fixing ui and taking expectations with respect to B. 
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3 Estimates on the solutions 

In this section, we derive moment estimates for the solutions to equations of the form: 



d 

dXt = VoiXt) dt + ^V,{Xt)odBi , t>0, (3.1) 



where x = Xo e M" and Vi,i ^ 1,2, ■■ ■ ,d are bounded, C°° vector fields in M" with 
bounded derivatives, and Vq is a possibly unbounded C°° vector field with bounded 
derivatives. Our bounds will be purely pathwise bounds, so the fact that the i3"s are 
sample paths of fractional Brownian motion is irrelevant, except to get a bound on their 
Holder regularity. 

AU that we will assume is that the driving process is 7-Holder continuous for some 
7 > 1/2. The main result in this section are a priori bounds not only on (3.1), but also 
on its Jacobian and its second variation with respect to the initial condition. This will 
be a slight generalisation of the results in [HN07], which required the drift term Vq to 
be bounded. 

Let ipt{x, w) = Xt denote the flow map solving (3.1) with initial condition Xq = x 
and define its Jacobian by 

_ dLpt{x,uj) 

■Jo,t = 5 ■ 

ox 

For notational convenience, set V = (Vi, V2, ■ ■ ■ , Vd). Then, we can write (3.1) in 
compact form: 

dXt = Vo(Xt) dt + V(Xt) o dBt , (3.2a) 

dJo.t = V^{Xt)Jo,tdt + V'{Xt)Jo.t o dBt , (3.2b) 

dJ^} - -J-lV^(Xt) dt - JoJV'(Xt) o dBt . (3.2c) 

Here, both J and J^^ are nx n matrices, and Jo,o = Jqo = 1. One crucial ingredient 
in order to obtain the bounds on the Malliavin matrix required to show the strong Feller 
property is control on the moments of both the solution and its Jacobian. These bounds 
do not quite follow from standard results since most of them require that Vq is also 
bounded. However we cannot assume this since we need condition (1.2) for showing 
ergodicity. 

3.1 Moments estimates 

Let CkiXi, . . . ,Xi^;y) denote the set of fc-multilinear maps from Xi x . . . x X/^ to 
y. As usual, Ci is denoted by C. In this section, we consider a system of differential 
equations of the form 

xt=xn+ Mxr)dr+ f(xr)dBr, (3.3a) 



Jo 
t 



yt^yo+ / gixr)(yr, dBr) , (3.3b) 







where xt G M", yt S R™, B : [0, co) i-^ M'' is a Holder function of order 7 > 1/2 
and /o: M" W\ f : ^ C{W^-W% g : W ^ £2(1^™, R''; R") are given 
functions. Note that (3.2) is indeed of this form with m = Ir?. 
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Remark 3.1 The reason why we do not include a 'dr' term in (3.3b) is because this 
is already covered by the present result by setting Bi{r) = r for example. Treating 
this term separately like in (3.3a) might allow to slightly improve our results, but the 
present formulation is sufficient for our purpose. 

Regularity of solutions to equations of this kind has been well studied, pioneered 
by the technique of Young ([You36]) and more recently using the theory of rough 
paths (see for example [LQ02, FrilO] and the references therein). Using fractional 
derivatives, it was shown in [HN07, Thm 3] that 

sup \yt\<2'+''^\<(^r,\y,\, (3.4) 

0<t<T 

where |]i?||o,T,7 is the Holder norm defined in (2.3) and M is a constant depending 
on the supremum norms of /q, /, g and their first derivatives. As mentioned earlier, 
these estimates are not sufficient to obtain moment bounds on the Jacobian Jg.t and 
its inverse, since we wish to consider situations where /q is unbounded, so that its 
supremum norm is not finite. 

However it turns out that an estimate of the type (3.4) holds even if only the deriva- 
tive of /o is bounded, thanks to the fact that the corresponding "driving noise" (t) is 
actually a differentiable function. This is the content of the following result: 



Lemma 3.2 For the processes Xt and yt defined in (3.3), we have the pathwise bounds: 



<A^(l + |a:o|)(l + ||S||o,T,7) 



1/7 



||?/||o,T,7<^/(l + l^-o|)|yo|e*'"""'''--, 



(3.5a) 
(3.5b) 



where M is an absolute constant which depends only on WfWao, ||/'||oo> ||/ol|oo> ||.9||oo> 
1 1 (7' I loo and T. 

Proof. The proof is almost identical to that of [HN07], so we only try to highlight the 
main differences. For s, t € [0, T], we have the apriori bound (see [HN07, p. 403]), 



f{Xr) dBr 



< M\\B\\o^T^^{(t - s)T + \\x\U,tMt - s)2t) 



where the constant M is indepedent of |li?||o,T,7, but depends on /. Since /o grows at 
most linearly, we also have 



faixr) dr 



< M{{t ~s)+ \xs\{t - s) + ||x|U,t,7(t - s)^+i) , 



and therefore 



sup 

s<m,n<t 



Thus we have 



foixr) dr 



< M 



((t -s)+ \xs\{t -s) + 2\\x\\s,tMt - s)''-^') ■ (3-6) 



\x\\s,tn < (1 + l|S||o,r,7 + - + (PIIo,t,7 + DiklU.t.iC* " ^V) 
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for a constant AIi depending on / and the terminal time T. Set A = (2Mi(||i3||o,T,7 + 
1))"^^'^. Then for |t-s| < A, 

\\x\\sj,y < 2A/i(l + ||B||o,T,7 + kslAi--^) (3 7^ 

Since by definition we have: 

||a;|U,t,oo < l^^sl + ||a;||s,t,-y|i - s|'' , 
for |t - s| < A, from (3.7) it follows that 

||a^|lM,oo < ks|(l + 2MiA) + l . (3.8) 
Iterating the above estimate for N ^ T/A, it follows that 

N-l 

Ikllo.T.oo < l^old + 2AfiA)^ + 1^(1 + 2MiA)'= < (1 + |xo|)(l + 2AfiA)^7V 

k=l 

= (l + |xo|)(l + 2AfiA)^/^r/A. 
Using the fact that (1 + ^)^ < for A > 0, this finally yields 

||^c||o,T,co < (1 + \xo\)e^^''^T(2Mi{\\B\\o.T.^ + 1)^'^ . (3.9) 

Note that in the last step of the above argument we bounded (1 + 2MiA)^/^ 
from above by e^^^^-'" and thus obtained a constant which is exponential in T, but 
independent of A = OdlBUg y''!^). If the noise corresponding to the vector field /o 
were not differentiable but only Holder continuous, then instead of (3.8), the exponent 
would include a power of the Holder constant of the driving noise, thus yielding an 
estimate comparable to those obtained in [HN07]. 

Thus substituting the bound (3.9) in (3.7) yields that, for t — s < A, 

\\x\\s,t,^ < M{1 + |a;o|)(l + ||B||o.t.7) • (3.10) 

Using the fact the function f{x) = x'' is concave for 7 < 1, it is shown in Lemma A. 2 
(proof given in the Appendix) that for = < mi < M2 < • • ■ < mat-i < un ^ T, 
we have the bound 

||a;||o.T.7 < max ||a;|U,.„,+,,^ . (3.11) 

Thus from (3.10) and (3.11) we deduce that 

lk||o,T,7 < ^'^(1 + ko|)(l + ||B||o,T,7)^'^'' = Mil + |xo|)(l + \\B\\o.T.,j)'^'' , 

proving the claim made in (3.5a). 

Now we come to the second part of Lemma 3.2. Once again, for s, t e [0, T] we 
have the apriori bound (see [HN07, p. 406]) 

g(Xr)yrdBr 



< Mil + \\B\\o.Tn)[\\9\\oo\\y\\s,tMt - (3.12) 

+ i\\g\U\y\\s,tn + y\U\y\\s,t.oo\\x\\s,tn)(t - ^f") 
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for a constant M independent of B. Thus, setting as before A = {2Mi{\\B\\o.t.-i + 
1))"^^'^, we deduce from (3.7) and (3.12) that, 

\\y\\s,t,j < M2(l + ||Bl|o,T,7)(ll?/||s,f,co+(l+||S|lo^T.7)lly|lM,ooA"'+|l2;|U,t^^A^) , 

(3.13) 

for a large enough constant M2 > Af(l + \xo\Te^^'^'^ + \\g\\oo + ||.9'||cc)- From the 
estimate (3.13), a straightforward calculation (for instance, see [HN07, p. 407]) yields 

||y||o,T,oo<A/2™lt^ilo--)^^^|yo|, 
bllo.T,, < Mil + l|Sl|o,T,,)^/^2*^^(i+ll^ll°.--'^^^|yo| , 
for a large constant M, thus concluding the proof. □ 

In the next sections, we will also require bounds on the second derivative of the so- 
lution flow with respect to its initial condition. Given an equation of the type (3.3), the 
second variation zt G M*' of x with respect to its initial condition satisfies an equation 
of the type 

Zt = Zo+ / h{Xr){Zr, dB,.) + / h{Xr)iyr, Ur, dBr) ■ (3.14) 

Jo Jo 

Here, the maps h: W C2(W,R'^;W) and h: W -> £3(1^'", K'", K''; K^) are 
bounded with bounded derivatives. Bounds on z are covered by the following corollary 
to Lemma 3.2: 



CoroUary 3.3 For the process zt defined in (3.14), we have the pathwise bounds: 

\\z\\o,Tn < Mid + \xo\')a + I2/0I') + kold + \xo\)) exp(A/||B|lJ/^^^) , 

where M is an absolute constant which depends only on \\f Woo, ||/'||oo> ||/ol|oo> ||.9lloo> 
Wg'Woc andT. 

Proof. Denote by Zt the '^'-valued solution to the homogeneous equation 

Zti^^+ f h{Xr){Zr£., dBr) , i^W . 

Jo 

Note that the inverse matrix Z^T^ then satisfies a similar equation, namely 

^r'e = e - / Z-^h(Xr)it dBr) , ^ e MP . 
"'0 

It follows from Lemma 3.2 that both Z and Z^^ satisfy the bound (3.5b), just like y 
did. On the other hand, z can be solved by the variation of constants formula: 



Zt = ZfZo + Zt Zg ^h{xs)iys, Us, dBs) 



(3.15) 



We now use the fact that if / and g are two functions that are 7-Holder continuous for 
some 7 > |, then there exists a constant M such that 



mdg(t) 



0,T,7 



<M(||/||o,T,7 + 1/(0)1) ||5||o,T,7- 



(3.16) 
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(See [You36] orfor example [MN08, Eqn2.1] for a formulation that is closer to (3. 16).) 
Furthermore, 

ll/.9llo,T,7 < Af(|/(0)| + ||/!|o,T,7)(l<?(0)| + ||.g||o.T,^) , 
so that, from (3.15), 

+ C(1 + ||Z|| )(1 + 11^-^11 O.T.-y) 

X ||/i'||oo(|a;o| + ||a;||o,T,7)(|j/o| + ||y|lo,T,7)^ll^llo,T,7 ■ 

The claim now follows by a simple application of Lemma 3.2. □ 

The above two lemmas immediately give us the required conditional moment esti- 
mates for the Jacobian Jo,t- 

Lemma 3.4 Let Xt, Jq t,mt be as defined in (3.2a), (3.2b) and (2.8) respectively. 
Then for any T > 0, p > 1 and 7 € (1/2, H) we have 

]E(ll>/o,.|IS,T,7) < A'^e*^*'+'l""°'--)''\ 

lE(l|Jo:'llo,T,7) <^e*^^'^"""°'""^''' ■ 
(Recall that E was defined in Remark 2.8.) 

Proof. We first write Bt ~ Bt + nit as in (2.8). Thus (3.2b) can be written as: 

dJo,i = V^(Xt)Ja.tdt + V'(Xt)Jo,t o (dBt + dmt) . 

By applying Lemma 3.2 to the SDE in (3.2b) with with yt = (t, Bt + rrit), fo = Vq, 
f ^ (0,Vi, . . . , Vd), g = (Vq, yj) and using the fact that I/7 > 1 we obtain: 

^ Me*^(i+ll^llo'^,.+ll™li;f^..) . (3.17) 

The fact that 7 > ^ is crucial. Since i? is a Gaussian process, it follows by Fernique's 
theorem [Bog98] that the Holder norm of B has Gaussian tails, thus proving our first 
claim. An identical argument for Jq^ finishes the proof. □ 

Remark 3.5 In an identical way, one can obtain similar bounds on the conditional 
moments of the solution and of its second variation with respect to its initial condition. 

4 Smoothness of conditional densities 

In this section, we consider again our main object of study, namely the SDE 

d 

dXt = VdXt)dt + Wt) o dBi , = .To e R" , 

1=1 

where Bt is a two-sided fBM. Our aim is to show that, if the vector fields V satisfy 
Hormander's celebrated Lie bracket condition (see below), then the conditional distri- 
bution of Xt given {Bg : s < 0} almost surely has a smooth density with respect to 
Lebesgue measure on M" for every i > 0. 
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In [BH07], the authors studied the same equation and showed that the law of Xt has 
a smooth density, but with the key difference that the driving noise Bt was a one-sided 
fBm and no conditioning was involved. Being able to treat the conditioned process will 
be crucial for the proof of the strong Feller property later on, this is why we consider it 
here. As explained in (2.8), we can rewrite Bt as Bt = Bt + rrit and thus obtain 

dXt = Vo(Xt) dt + V(Xt) dmt + V(Xt) o dBt . (4. 1) 

While the distribution of the conditioned process Bf is different from that of the one- 
sided fBm, its covariance satisfies very similar bounds so that the results from [BH07] 
would apply. (See Proposition 4.4 below.) 

However, the difficulty comes from the time derivative of the conditional mean mj, 
which diverges at t = like 0{t^~^). To use the results from [BH07] directly, we 
would need the boundedness of the driving vector fields, uniformly in time, so that 
they are not applicable to our case. Our treatment of this singularity is inspired by 
[HM09], where a related situation appears due to the infinite-dimensionality of the 
system considered. Using an iterative argument, we will show the invertibility of a 
slightly modified reduced Malliavin matrix which will then, by standard arguments, 
yield the smoothness of laws of the SDE given in (4.1). 

Before we state Hormander's condition, let us introduce a shorthand notation. For 
any multiindex I — {ii,i2, ■ ■ ■ , ik) G {0, 1, 2, • • . , d}'^, we define the commutator: 

Vi = [v,,AV^,,--- AV^,_,,v,J ■■■]]. 

We also define recursively the families of vector fields 

n 

Vn = {Vi,I e{0,--- ,dr-^ X {!,■■■ ,d}} , K^IJVfc. (4.2) 

fe=i 

With these notations at hand, we now formulate Hormander's bracket condition [Hor67] : 

Assumption 4.1 For every xq S M", there exists N £ N such that the identity 

span{U{xo) : U € Vat} = R" (4.3) 

holds. 

4.1 Invertibility of the Malliavin Matrix 

It is well-known from the works of Malliavin, Bismut, Kusuoka, Stroock and others 
[Mal78, BisSl, KS84, KS85, KS87, Nor86, Mal97, Nua06] that one way of proving 
the smoothness of the law of Xt (at least when driven by Brownian motion) is to first 
show the invertibility of the 'reduced Malliavin matrix' ' : 

Ct = AtA*t^ / JolV{X,)V{X,riJolT ds , (4.4) 
Jo 

where we defined the (random) operator: At ■ ^^([0, T], M'^) i-^ M" by 

Atv= [ JolV(Xs)v(s)ds . (4.5) 

' This is a slight misnomer since our SDE is driven by fractional Brownian motion, rather than Brownian 
motion. One can actually rewrite the solution as a function of an underlying Brownian motion by making 
use of the representation 2.8. but the associated Malliavin matrix has a slightly more complicated relation to 
C't than usual. This will be done in Theorem 4.1 1 below. 
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Unfortunately, for the stochastic control argument used in Section 5 below for the proof 
of the strong Feller property, the invertibility of Ct turns out not to be sufficient. We 
now define a family of smooth functions hx ■ [0, T] -> [0, 1] : 

, , , f 1 if s < T/2, 

^^(^^=1 if,s>3T/4. ^^-^^ 

and consider the (random) matrix Ct given by 

Ct = At HtAt ■ (4.7) 

(See Remark 5.5 below for the reason behind introducing the function Ht and the 
matrix Ct ) Obviously, Ct < Ct as positive definite matrices, so that invertibility of 
Ct will in particular imply inertibility of Ct, but the converse is of course not true in 
general. Define the matrix norm: 

\\Ct\\ = sup {v,Ctv), ipeU" . 
Since Ct is a symmetric matrix, by the spectral theorem we have 

WCt'W^—. — ^v^' 

inf|<p|=i(u, Ctw) 

The aim of this section is to show that Ct is invertible, and to obtain moment estimates 
for II C^^. The main result of this section can be formulated in the following way: 

Theorem 4.2 Assutne that the vector fields {Vi \ have bounded derivatives of all orders 
for i > Q and are bounded for z > 0. Assume furthermore that Assumption 4.1 holds, 
fixT > and let Ct be as in (4.7). Then, for any R > and any p > 1, there exists a 
constant M such that the bound 



P inf {ip, Ct^) <e] < Me^ , p &W , 
\||<^||=i / 

holds for all e G (0, 1] and for all xq and uj such that |a;o| + ||cj||(^^a) < R. Here, xq is 
the initial condition of the SDE (4.1). 

Corollary 4.3 The matrix Ct is almost surely invertible and, for any R > and any 
p > I, there exists a constant M such that 

EllCy^ll <M , 

uniformly overall xq and uj such that \xq\ + ||a;||(^ ,5) < R. 

Proof. Notice that for p > 1, we have HCj^'^H = HC^^p since Ct is symmetric an 
positive definite. Thus 

/•oo 

E\\Cj,P\\ =p xP-^FiWCj^^W > x)dx 



p ( xP-ipf inf (p, V) <-\dx<oo, 

Jo V||v>j|=l X/ 



where the last inequality follows from Theorem 4.2, since for a; > 1 there exists a 
constant (depending only on i?) such that P( inf||^ 11^1 ((y9, Cj^^(p) < j) < Mx~^~p 
and we are done. □ 
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Before proceeding further, we need the following version of Norris' lemma [KS85, 
Nor86] adapted from [BH07] which will be the key technical estimate needed for show- 
ing the invertibility of Ct- In the case of white driving noise, Norris' lemma essentially 
states that if a semimartingale is small and if one has some a priori bounds on its com- 
ponents, then both the bounded variation part and the martingale part have to be small. 
In this sense, it can be viewed as a quantitative version of the Doob-Meyer semimartin- 
gale decomposition theorem. One version of this result for fractional driving noise is 
given by: 

Proposition 4.4 Let H G (5,1) and 7 g and let B be as in (2.8). Then there 

exist exponents q > and rj > such that the following statement holds. 

Consider any two families a^, 6^ of processes taking values in M and R'' respec- 
tively such that V,{\\a'^\\^ j, ^ +\\b'^\\^ j, ^) < Afp£"''P, uniformly in e € iQ^X], for some 
Mp > and for every p> 1. Let be the process defined by 

yt = f alds+ f blodBs, te[0,T] . 
Jo Jo 

Then, the estimate 

< £ and |la||o,T,oo + |l&|lo,T,oo > e"^) < C^e^ 
holds for every p > 0, uniformly for e G (0, 1]. 
Proof Setting f(s, t) = E\B(t) - B(s)\^ as in [BH07], we have 

(We use the convention s < t.) A somewhat lengthy but straightforward calculation 
shows that this can be rewritten as 

/(., t) = It - + ((1 + - dx) . 

Similarly, one has 

i-s/\t-s\ 

Jo 

Since the integrands of both of the integrals appearing in these expressions decay like 
^2H-3^ they are integrable so that B is indeed a Gaussian process of 'type H\ see 
[BH07, Eqn3.5]. 

The proof is now a slight modification of that of [BH07, Prop 3.4]. It follows from 
[BH07, Equ 3.18] and the last displayed equation in the proof of [BH07, Prop 3.4] that 
there exist exponents q and f3 such that the bound 

P{\\y\\o,T,oo < e and ||a||o,T,co + ||&||o,t,co > e"^) < Cpe" + P(||6||2 + ||a||2 > s'^) 



holds. The claim now follows from a simple application of Chebychev's inequality, 
provided that we ensure that ij is such that 27] < p. □ 



Smoothness of conditional densities 



16 



For making further calculations easier, we now introduce some definitions that are 
inspired by the proofs in [HM09] . 

Definition 4.5 A family of events A = {j4£}£<i of a probability space is called a 
"family of negligible events" if, for every p > I there exists a constant Cp such that 
IP(^e) < CpSP for every e < 1. 

Definition 4.6 Given a statement $e, we will say that " $e holds modulo negligible 
events", if there exists a family {A^} of negligible events such that, for every e < 1, 
the statement $e holds on the complement of A^. 

With the above definitions. Theorem 4.2 can be restated as: 
Tlieorem 4.7 Fix R > 0. Then, modulo negligible events, the bound 

inf (os, Ct'-p) > e 
llvll=i 

holds for all xq and lo such that \xa\ + \\Lo\\{y_s) < -R- 

Remark 4.8 The family of events in question can (and will in general) depend on R, 
but once R is fixed the bounds are uniform over [xol + .5) < R. While this 
uniformity is not really required in this section, it will be crucial in Section 5. 

Proof of Theorem 4.7. It suffices to show that for any fixed ip with \\(p\\ = 1, the bound 
{(fi, Ct'p) > £ holds modulo a family of negligible events. The conclusion then fol- 
lows by a standard covering argument, see for example [Nor86, p. 127]. Furthermore, 
throughout all of this proof, we will denote by M a constant that can change from one 
expression to the next and only depends on R and on the vector fields V defining the 
problem. 

Fix now such a (p and assume "by contradiction" that {ip, Crf) < £■ We will show 
that this assumption leads to the conclusion that 

sup \{'f,U{xo))\<e'^ , for some g'>0, 

modulo negligible events. However, if Hormander's condition is satisfied, the conclu- 
sion above cannot hold since {U{xq), U G Vat} span M" so that there exists a constant 
c such that supj^^y^ | (cp, U{xo)) \ > c almost surely, thus leading to a contradiction 
and concluding the proof. For any smooth vector field U, let us introduce the process 
Zu(t) = J^ jUiXt). With this notation, we have 

Lemma 4.9 There exists r > such that, modulo negligible events, the implication 

{(p,CT(p)<e =^ i|(v5,^c/(-)>llo,T/2,oo < 
holds for all U £ Vi. 

Proof. By the definition of h, we have the bound 

{<f, Ctv) = V / {if, J^lV,(X,)f hris) > V / J^]V,{X,)fds . 

(4.8) 
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We now wish to turn this into a lower bound of order (for some r > 0) on the 
supremum of the integrand. Our main tool for this is the bound 

||/||o,T,oo<2max(r-i/2||/j|^,j„,,j, II/IIII^tj H/HoT;) ' 

which holds for every 7-Holder continuous function / : [0, T] n- M. The proof of this 
is given in the appendix as Lemma A. 3. 

By Lemma 3.4 and Chebyshev's inequality it follows that || {ip, Jq^V(X.)) \\o,t,i < 
£^^/^ modulo negligible events. Thus if {ip, Ct'-p) < £, from Lemma A.3 and (4.8) we 
deduce that, for each 1 < j < d 

\\{^, Jo'V,(X.))\\o,T/2,o. < Af£^ , (4.9) 
modulo negligible events. The claim follows by choosing r < • '-' 

We will now use an induction argument to show that a similar bound holds with Vj 
replaced by any vector field obtained by taking finitely many Lie brackets between the 
Vj's. Whilst the gist of the argument is standard, see [KS85, Nor86, BH07], a technical 
difficulty arises from the divergence of the derivative of m(t) at t = 0. To surmount this 
difficulty, we take inspiration from [HM09] and introduce a (small) parameter a > 
to be specified at the end of the iteration. Then, by Lemma 2.5, we have 

^ < Af£"<^-i' . (4.10) 

at e»,T/2,oo 

Instead of considering supremums over the time interval [0, T/2], we will now consider 
supremums over [e°',T/2] instead, thus avoiding the singularity at zero. 
In order to formalise this, for ?- > and i > 1, define the events 

Kl(a,r) = {\\{ip, Zu(Me^,T/2,o. < e\ V C/ e V,} . 

With this notation at hand, the recursion step in our argument can be formulated as 
follows: 

Lemma 4.10 There exists q > such that, for every a > and every i > 1, the 
implication 

K^(a, r) ^ K^+i (a, gr + a{H - 1)) , 
holds modulo negligible events. 

Proof. Assume that K|(q!, r) holds. For any t < T/2, it then follows from the chain 
rule that 

{ip,Zu(t) - Zu(e")) = Jjifi, Jo^llVo, U](Xs)) ds (4.11) 
+ J2 [ Jo,li^^, U]{X,)} dml +J2 [ iv, Jo,l^V^^ U]{Xs)) o dBl , 

i=l i=l 

where 13. is as in (2.8). 

Now by the hypothesis. Proposition 4.4 and (4. 1 1), it follows that for any U G Vi, 



d 



< Afe'"9 , (4.12a) 

e°,T/2,oo 
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J^hVo, U](X.) + ^JohV^, U](X.) 

i=l 



< Me'''' , (4.12b) 

e°,T/2,oo 



modulo negligible events, with q as given by Proposition 4.4. Plugging (4.10) and 
(4.12a) into (4.12b), it follows that for every U G Vi, 

\\{^,Jo'[Vo,U](X.))\\e-,T/2,oo < Me'-'^+"(^-i» 

modulo negligible events. Choosing any q < q, the claim then follows from the fact 
that for every exponent S > 0, one has Me^ < 1 modulo negligible events. (Recall 
that < 1, so that (4. 12a) is actually slightly better than the claimed bound.) □ 

Now we iterate the previous argument. Let iV G N be such that 

inf inf V \{^M(x))\^ > . (4.13) 
bl<«l?l=i 

It follows from the Hormander condition and the fact that the vector fields are smooth 
that such an N exists. From Lemmas 4.9 & 4.10, we obtain that for every U € Vjv, 
modulo negligible events, 

||(</7, Jo-.i[/(X))||,„,T/2,oo < e'^ , (4.14) 

where qi = rq^ + a(H — l) ''^7-r^ ^ - Since q and r are fixed, we can choose a small 
enough such that qi > 0. 

Now we make use of the Holder continuity of Jq^ and X. Indeed, it follows from 
Lemma 3.4 that 

\\Jo'U(X.) - f/(xo)||o,s°,oo < , (4.15) 

modulo neghgible events, for any q2 < a-f. Combining (4.14) and (4.15), we conclude 
that 

sup \{ip,U{xo)) \ < e"' +6"" , 

U&Vn 

modulo negligible events. This is in direct contradiction with (4.13), thus concluding 
the proof of Theorem 4.7. □ 

We conclude that the matrix Ct defined in (4.4) is almost surely invertible and 
£(11(7^11"^) < M forp > 1. As a consequence of Theorem 4.2, we obtain the smooth- 
ness of the laws of Xt conditional on a given past uj G W: 

Theorem 4.11 Let ftix, uj) be the solution to (1.1) with notations as in Section 2.1, 
and assume that Assumption 4.1 holds. Then, for every u G W, the distribution of 
ift{x, cD) under V(id, ■ ) has a smooth density with respect to Lebesgue measure. 

Proof. In order to obtain the smoothness of the densities, it suffices by [Nua06, Mal97, 
BH07, NS09] to show that the Malliavin matrix Mr of the map W ^ X^, where 
W is the standard Wiener process describing the solution as in (2.8), is invertible with 
negative moments of all orders. As we will see in Section 5 below (see equation 5.3), 
A^T is given by 



L2 
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where the L^-adjoint of the fractional integral operator 2 is given by 

We already know from Theorem 4.2 that ||.4* Jq -^^1112 > £ modulo negligible events 
and we would like to make use of this information in order to obtain a lower bound on 
Mt- 

Since 2 and 2?^^ 2 are inverses of each other, for any measurable / we have, 

ll/lli. = {f,i"---v"-if)L. < ||(i^-3)*/||^.||2?«-^/!|i2 
<M||(i^'^)7IL2||/||o,T,,. 

Thus we obtain, 

|l(I^-^)*^VoVe!L. > M\\A\Jlr^\\h\\A*Jlrao.T,, ■ 

Since we know from the results in Section 3 that i|.4* Jo,tCI1o,t,7 < modulo 
negligible events, we conclude that 

\\il"-irA*J^,T^h.>e\\A*Jlrah, 
holds modulo negligible events. The proof now follows from Theorem 4.2. □ 

5 Strong Feller Property 

Recall from the arguments in Section 2.2 that the strong Feller property (2.2) is the 
only ingredient required for the proof of Theorem 1.3. We will actually show some- 
thing slightly stronger than just (2.2). Fix some arbitrary final time T > 1, a Frechet 
differentiable bounded map tp: C([l,r],R") — )- R with bounded derivative^, denote 
by Ri^T ■■ C(R+, M") C([l, T], M") the restriction operator, and set as before 

Qtp(x,uj)= / tlj{Ri^Tz)Q{x,uj;dz) . 

JC(R+,K<*) 

With these notations, the main result of this article is: 

Theorem 5.1 Assume that the vector fields {Vi\i>Q have bounded derivatives of all or- 
ders and that furthermore the vector fields {Vi\i>i are bounded. Then, ifHonnander's 
bracket condition (Assumption 4.1) holds, there exists a continuous function K : R" — >■ 
M_l- and a constant M > such that for lo G W^.a and x e R", the bound 

\D,Qij(x,io)\ < if(x)e^^"^'""of;.-, , (5.1) 

holds unifonnly over all test functions tp OS above with ||V'||oo < 1- Here, Q was defined 
in (2.6) and both K and M are independent of T andofip. In particular, (1.1) is strong 
Feller. 

^The space C([l, T], R") is endowed with the usual supremum norm 
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The remainder of this section is devoted to the proof of this resuh. For conciseness, 
from now onwards we will use the shorthand notation 

mx,uj)\<E(x,u;), 

to indicate that there exist a continuous function K on R" and a constant M such that 
the expression '^{x, uj) is bounded by the right hand side of (5.1). 

To obtain the bound for DQip(-,uj), we use a 'stochastic control' argument which 
is inspired by [HM06, HM09] and is similar in spirit to the arguments based on the 
Bismut-Elworthy-Li formula [EL94, DPEZ95]. We will explain this argument in more 
detail in Section 5.1 below, but before that let us introduce the elements of Malliavin 
calculus required for the construction. 

Let and Jqj be the solution and its Jacobian, as defined in equations (3.2a) and 
(3.2b) respectively with initial condition Xq ~ x. Let (fiti^, uj) ~ Xf denote the flow 
map as in Section 2.1. Henceforth we will use iptix^ Cj) and Xf interchangeably, with 
a preference for the first notation when we want to emphasise the dependence of the 
solution either on the realisation of the noise or on its initial condition. 

By the variation of constants formula, the Frechet derivative of the flow map with 
respect to the driving noise Cj in the direction of v{s) ds is then given by 

D^,X^ = D.„Lpt{x, Co) = Jo^AtV , (5.2) 

where At is as defined in (4.5). 

Remark 5.2 Notice that the operator D defined above is a standard Frechet deriva- 
tive unhke in the case H = ^ where it would have to be interpreted as a Malliavin 
derivative. This is because the integration with respect to the fBM for H > ^ is just 
the Reimann-Stieltjes integral, whereas for H = ^ the stochastic integral is not even 
a continuous (let alone differentiable) map of the noise in the case of multipUcative 
noise. 

As already mentioned in Remark 2.4, the operator 'D^~ ^ and its inverse i provide 
an isometry between the Cameron-Martin space of the Gaussian measure Viuj, ■ ) and 
that of Wiener measure. 

As a consequence, if v is such that v = V^^^y e L^(]R_|-, M''), then we have the 
identity 

D^Xf = D^iptix.Cj) = \tAtv = Jo,tAtl"-h , (5.3) 

where D is the Malliavin derivative of when interpreted as a function on Wiener 
space via the representation (2.8). 

Before we continue, let us make a short digression on the definition of the Malli- 
avin derivative that appears in the above expression. In general, if Z : >V — >^ M is any 
Frechet differentiable function, then its Malliavin derivative with respect to the under- 
lying Wiener process is given by 

n,Z = D,Z , (5.4) 

where D^Z denotes the Frechet derivative of Z in the direction v{t) dt and where v 
and V are related by u = V^^^v as before. 

Since, for every T > 0, the map T)^^ 2 restricted to functions on [0, T] is bounded 
from W to C([0, T], R'') and since the dual of the latter space consists of finite signed 
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measures, we conclude from (5.4) that if Z is Frechet differentiable, then there exists 
a function s M' DgZ which is locally of bounded variation and such that the Malliavin 
derivative of Z can be represented as 

/"OO 

DtZ= {D,Z,i){s))ds . 
Jo 

(The scalar product is that of here.) The quantity DgZ should be interpreted as 
the variation of Z with respect to the infinitesimal increment of the underlying Wiener 
process W at time s. 

With this convention, it follows from (5.3) and the definition of 2 that for some 
constant ch one has the identity 

D,Xf = CH f Jr,s V{X^) {r - t)"~i dr , (5.5) 

J s 

for s < t. Here, V denotes the matrix (Vi, . . . , Vd) as before. Note also thatDgXf = 
for s > t since X is an adapted process. Of course, the whole point of Malliavin 
calculus is to also be able to deal with functions that are not Frechet differentiable, but 
for our purposes the framework discussed here will suffice. 

5.1 Stochastic control argument 

At the intuitive level, the central idea of the stochastic control argument is as follows. 
In a nutshell, the strong Feller property states that for two nearby initial conditions .tq 
and j/o and an arbitrary realisation to of the past of the driving noise, it is possible to 
construct a coupling between the solutions xt and yt such that with high probability 
(tending to 1 as y x), one has xt = yt for t > 1. One way of achieving such a 
coupling is to perform a change of measure on the driving process for one of the two 
solutions (say y) in such a way, that the noise ojy driving y is related to the noise Ux 
driving x by 

didy(t) = dujxit) + Vx,y{t) dt , 

where v is a 'control' process that aims at steering the solution y towards the solution 
X. If one takes y = a; + and looks for controls of the form v^.y = ev, then in the 
limit £ 0, the effect of the control v onto the solution x after time t is given by (5.3), 
while the effect of shifting the initial condition by ^ is given by Jo,tC- It is therefore not 
surprising that in order to prove the strong Feller property, our main task will be to find 
a control v such that Jq iAiv = Jo.i^, or equivalently Aiv ~ ^. This is the content of 
the following proposition: 

Proposition 5.3 Assume that, for every x G K" and ui G W, there exists a stochastic 
processes v G C'''(M+, W^) such that the identity 

Av = i, 

holds almost surely, where A == Ai defined in (4.5). Assume furthermore that the map 
w 1-^ V is Frechet differentiable from W to C'(]R+, M''). Finally, we assume that the 
Malliavin derivative of the process v = 'D^~^v satisfies the bounds 



/•oo /"OO /"OO 

Ml = E / |w(s)p ds < 00 , M2 = E / / |Dtw(s)p dsdt < 00 . (5.6) 
Jo Jo Jo 
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Tlien, the bound 



holds uniformly over all test functions ijj as in Theorem 5.1, all T > I, and all ^ S M". 

In particular, if v is such that Mi + M2 < E(a;, ui), then the conclusions of Theo- 
rem 5.1 hold. 

Remark 5.4 Since V"- 5 is bounded from , W^) into C(M+ , R''), the MalUavin 

derivative of v exists, so that the expression in (5.6) makes sense. 

Proof. Given an initial displacement ^ e R", we seek for a 'control' v on the time 
interval [0, 1] that solves the equation Av = ^, where we use the shorthand notation 
^ = ^1 . If this can be achieved, then we extend v to all of by setting v{s) = for 
s > 1 and we define v = V^^^y. Note that since v{s) = for s > 1, it follows from 
the definition of At that we have Atv = Av = ^ for t > 1. If w is sufficiently regular 
in time so that this definition makes sense, we thus have the identity 

D^Xf; = Jo,tC , (5.7) 

for every T > 1. In particular, for a test function ip as above and for ^ G M", 

{D.Q^ix, w), = -E((Dij)(cp(x, lD)) Jo..c) 



It then follows from the integration by parts formula of Malliavin calculus [Nua06] that 
this quantity is equal to 



= ]E(v(¥'T(a;,c^)) / {v{s),dW,)^ < U'\ 



T 



The stochastic integral appearing in this expression is in general not an Ito integral, 
but a Skorohod integral, since its integrand is in general not adapted to the filtration 
generated by W. 

For Skorohod integrals one has the following extension of Ito's isometry [Nua06]: 

{i>{s),dWs)^ |f>(s)pds^ +E / / tT{Dtv{sfDsV{t))dsdt 

< Ml + M2 , (5.8) 

where as before D denotes Malliavin derivative with respect to the underlying Wiener 
process. The claim then follows at once. □ 

In order to have the identity (5.1), it therefore remains to find a control v satisfying 
Av = ^ and such that 

A/i +M2 < E(x,a;), (5.9) 

where the quantities AIi and A/2 are defined by (5.6) in terms of v — V^^^v. A nat- 
ural candidate for such a control v can be identified by a simple least squares formula. 
Indeed, note that the adjoint A* -.W^ ^ L'^(^+,W^) of A is given by: 

{A*Ois)^V{X,nj^lr£,, CeK"- (5.10) 
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Let h = hi where hi : [0, 1] 1R+ is as defined in (4.6). Assuming that AhA* is 
invertible, one possible solution to the equation Av = is given by 

v{s) EE h(s) {A*(AhA*r^O(s) ■ (5.11) 

Note that AhA* is precisely equal to the 'modified Malliavin matrix' Ci defined in 
(4.7). By the results from the previous section the matrix Ci is almost surely invertible 
and therefore the expression (5.11) does make sense. 

The aim of the next section is to show that the choice (5.11) does indeed satisfy the 
assumptions of Proposition 5.3. 

Remark 5.5 For the identity (5.3) to hold, the stochastic control v{-) needs to belong 
to the Cameron-Martin space of 7^(0, •), and the function h{-) is introduced just for this 
purpose. Indeed, since we extend the control v{-) to vanish outside the interval [0, 1], 
one needs to be careful about the regularity of v{s) at s = 1. It can be shown using the 
fractional integrals that, in order for v(-) to have the required regularity, the function 
h{s) should be ©((1 — s)^) for s sa 1. The specific /ii(-) we chose in (4.6), of course, 
satisfies this. 

5.2 Proof of the main result, Theorem 5.1 

In view of Proposition 5.3, the proof of Theorem 5.1 is complete if we can show that 
the stochastic process v defined by (5.11) satisfies the bounds (5.9). This will be the 
content of Propositions 5.8 and 5.11 below. 

We start with an estimate on the Holder norm of v which will be used repeatedly. 
Here and in the rest of this section, we will make use of the notation 

\M-I = ll«ilo,l,7 + IK'llo,l,oo 

for the Holder norm of v. This has the advantage of being a norm rather than just a 
seminorm and, for any two Holder continuous functions u and v, we have the bound 

for some fixed constant M > 0. With this notation, our bound is: 

Lemma 5.6 The stochastic control v defined hi (5.11) is continuous and satisfies the 
bound 

Ellwll^ < E{x,uj) . 
Proof. Setting v{s) ~ h(s)u(s), we obtain from (5.11) the bound 

ll"ll' < (5.12) 

where 77 = {AhA*)^^£,. From (5.10) and the fact that V has bounded derivatives it 
follows that 

\\A*v\U < M\\ViX.)\\JJ^^^\\^\r^\ < M||X||^|| Jo;.i||^|,7| , 
and therefore by Cauchy-Schwarz 
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From Lemma 3.4 we deduce that E(|lX|l^|lJo".^|l^) < E(.t,w). Also from Theorem 
4.2, we have E\(AhA*)^'^^\P < E(.t, lu) for any p > I. Thus we obtain 

]E||m||2 < E(x,u}) . 

The claim then follows immediately from v = hu. □ 

Remark 5.7 Had we not introduced the control h, we could still obtain continuity of 
on [0, 1], but we would have a discontinuity at t = 1. 

As a consequence, we show that v = V^^^v is square integrable: 

Proposition 5.8 For the stochastic control v(-) defined in (5.11) and the operator T> 
defined in (2.4), the following estimate holds: 

Before we turn to the proof of Proposition 5.8, we state the following useful result: 
Lemma 5.9 For any s > 0, a > and f <E with 7 > a , 

r(l - a) V^fis) = / [s- rr"f(r) dr 
as Jo 

= s-"f(3)-af (s-rr""i(/(r)-/(s))dr . (5.13) 
Jo 

Proof. See [SKM93, Thm 13.1], as well as Lemma 13.1 in the same monograph. □ 

Remark 5.10 The right hand side of equation (5.13) is called the Marchaud fractional 
derivative. Lemma 5.9 actually holds true for a larger class of functions /, namely 
those such that / = I" ,9 for some g € L^([0, 1], K). 

Proof of Proposition 5.8. The main idea behind this proof, and a recurrent theme in 
this section is the following. Although the stochastic control v{-) vanishes outside 
[0, 1], the integrand v{s) in non-zero for all s <S [0, 00). Thus the square integrability 
of v{s) on the interval [0, 00) is shown by dealing with the singularities of the integrand 
differently at s = 0, and s = 00. At s = 0, we use the fact that v is Holder continuous 
whereas at s = 00, we use the fact that v is bounded. 
By definition, 

d /" ^ 

v(s) = — (s — r)"^^^ v(r) dr . 



Applying Lemma 5.9 with a = £f — i, we obtain the bound 

\v{s)\<Ms^~"\v{s)\+M [ {s~r)-^-"\v{r)-v{s)\dr , 



for some constant M > 0. Noting that the control v(-) vanishes outside the interval 
[0, 1], we have 



E( / \iis)\''ds) < ME\\vf [ s^~^"ds (5.14) 
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ME f f {s-rr^-"\v{r)-v{s)\dryds 



The first integral s^~^^ ds in equation (5.14) is clearly finite since iJ < 1, so that 
this term is bounded by E(a;, uj) by Lemma 5.6. Now we show that the second term has 
a similar bound. To this end, set 

C(s)= [ (s-rr^-"\v(r)-v(s)\dr . 
Jo 

For s < 1, we obtain 

|C(s)H [ (s-rr^-"\v{r)-vis)\dr < M\\v\\^ [ g-^-^+'^ds 
Jo Jo 

< ]\I\\v\\^si~"+-^ , 

so that the required bound follows again from Lemma 5.6. For s > 1, we observe that 
C(s) has a square integrable singularity, and thus the boundedness of the control v(-) 
will be sufficient. Indeed it follows from the definition of h (4.6) that 

/•3/4 

IC(s)l < / (s~3/4:)-^-"\v(r)-v(s)\dr<M\\v\\^(s-3/4)-^-". 
Jo 

Since this expression is square integrable, it follows that K \(^{s)\ds < E{x,uj), 
and the claim follows from (5.14). □ 

We conclude with the corresponding bound on the Malliavin derivatives of the con- 
trol v: 

Proposition 5.11 Let v be as defined in (5.11), the operator T> as defined in (2.4) and 
D be the Malliavin derivative. Then following estimate holds: 

ti^j^ \\\)Ms)f dsdt^ <¥.{x,oj) . (5.15) 

Proof. The main difficulty here is to obtain control over the modulus of continuity of 
the quantity \)sv{t) for all s, t G [0, 1]. On route to establishing the Holder continuity 
of the Malliavin derivative applied to the stochastic control, we need to perform several 
careful estimates using the fractional derivatives. 

First notice that, since v{s) G T\ for every s > (J^t, t > denotes the filtration 
generated by the increments of W), we have Dtu(s)) = for i > 1, so that the inte- 
grand in (5.15) vanishes outside the set [0, 1] x [0, 1]. Furthermore, it is straightforward 
to check that Dt and T>^~ 2 commute, so that 



/* 00 /*00 /* 1 /* 1 

y \^tv{s)f dsdi)i <t[^j J \\{V"-iDtv){s)fdtds 
By Lemma 5.9, the claim follows if we can show that 

2 



'iff f is-r)--'-"(Dtvis)-Dtvir))dr 
^ Jo Jo Jo 



dtds] < E{x,uj) . (5.16) 
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Expanding the integrand, we obtain 

Dt ( /" (s - r)-i""{v{s) - v{r)) dr^ 



< I (s~ry^-"\Dtv{s)\dr 



+ 
< |s5 



H 



(s - ryi^" {v{s) - v{r)) drj 



{s-tr^-"\\\Dtv\\o,i.oc 
^ (s - rr^-"Dt{v(s) ~ v(r)) dr 



(5.17) 



Since the double integral of \s^^^ — {s — t)^^^\'^ is finite, the first term will be 
bounded by E(.t, uj) if we can show that E sup^^jQ ^ j|Dtu||o i oo < u). Since this 
bound will be a byproduct of the analysis of the second term, we postpone it to (5.23) 
below. 

The second term of equation (5.17) needs more delicate calculations. We will show 
that there exists a random variable K and an exponent a < 1 + j — H such that the 
following bounds hold: 

\T)tv(s)-Dtv(r)\<K\s-r\''\r~t\-°' , with E A'^ < E(a:, w) . (5.18) 

Assuming these bounds, we proceed similarly to the proof of Proposition 5.8: 

nl /.I 



E 



JQ 



is - rr^-"Dt{v{s) - vir)) dr 
■ 1 pi 



dtds 



< E K 



JO 
1 rl 



dt ds 



(s_r)-5-^+T(r-t)-"dr 

<e(^K^ (s-t)i+2''~^^~2"dtds) <Eix,uj), 

where we used the bound on a to obtain the last inequality. It therefore remains to 
show (5.18). 

The Malliavin derivative of v is given by 

T)Ms) = h{s)T)t(v(x,nj^inAhA*r^^) . 

Applying the product rule, we break this into a sum of three terms: 

DMs) = his) [-Dt {v{X,r){Jo.lnAhA*)-'^ + ViXsTDt {{Jo,lr)(AhA*r'^ 

+ ViXsTiJo^lTMUhA*)-'^) 

"li+Ia+Is- (5.19) 

To make the proof compact, we introduce the stochastic process given by col- 
lecting the various objects appearing in this expression: 



Y, = {Xs,V{X,)AJolT 



(5.20) 



It the follows from (1.1) and the expression (3.2b) for the Jacobian Jo t that the compo- 
nents of Ys solve equations of the form (3.5). One can check that both the Jacobian j7o,s 
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for Y and its inverse JTq I then satisfy equations of the type (3.14). So by Corollary 
3.3, we obtain 

Wh + WJoAU + \\Jo:'\\-, < K{xo)exp{M\\B + 0^1]]/-^) , (5.21) 

for some constant M and some continuous function K. 

Furthermore, just as in (5.5), the Malliavin derivative of Y is given by 

T>tYs Ja,s U{Ya) (a - tf-^'^ da lt<s , (5.22) 

for some smooth function U which grows linearly at infinity. 

We now tackle the bounds on the two very similar terms Ii, I2. Since t^^^ is 
integrable at and since U grows linearly at infinity, it follows immediately that 

sup |D,y,| < A/||Jo,.|Ml^o".'ll7ll^-il7 • (5-23) 

s,t6[0,l] 

Thus by (5.21) we have E supjgjQ l|Dti'|lo,i,cx) E(x,w). Regarding its modulus of 
continuity, note that 

DtY,~DtYr ^ Ja,sU(Ya)(a-t)"-Ua + {Ja,s-Ja,r)U(Ya)(a^t)''-i da . 

The first term is bounded by 

M\r~t\"-i\r~s\\\JoJJJo:%\\Y\\,, 
and the second term is bounded by 

M|r-sn|Jo,.|Ml^o:.'ll7ll>^ll7' 
so that there does exist a constant M such that 

\DtYs ~DtYr\ < M\r - t\"-l\r - s\\\JoA\M:X\\Y.h , 

uniformly over t < r < s < 1. Since 7 > i, we do indeed have | — i/< 1 + -f — H, 
so that this bound is of the type (5.18) with KK^ < E(x, w) for every p > 0. 

Notice furthemore that by Theorem 4.2, the matrix Ci = AhA* G M" is almost 
surely invertible and E|^/i^* j^^ < E(a;, to) for all p > 1. Combining this with (5.21) 
and the expressions for Ii and I2, we conclude that (5.18) does indeed hold for these 
terms. 

Let us now turn to the term I3. Recall that the Frechet derivative of the inverse an 
n X n matrix A in the direction _ff is given by D^^^^ = —A^^HA^^. It thus follows 
from the chain rule that 

Dt{iAhA*r^O = -(AhA*)-^ {Dt{AhA*)){AhA*r^^ 

Thus by Cauchy-Schwartz we have, 

(E\\mAhA*r'oC,^Jj' <nuhA*rr'^ 
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By the product rule and the definition of A, 



DtiAhA*) = / {OtJol) ViX,) his) V*{X,) (J^l)* ds + ... , (5.24) 
Jo 

where ". . ." stands for the remaining three terms that are very similar. Thus from (5.23) 
it follows that 



E 



{T)tJo,l)V(Xs) his) V*iX,) iJ^l)* ds 

< MEIlJo.ll^ llJoT-'ll^ ll^^ll'" \\ViX.)\\f < Eix.uj). 



Since J ^ and ViX) play indistinguishable roles in (5.20) and (5.21), the bounds on 
the remaining three terms follow in an identical fashion, so that 

E sup \DtiiAhA*)-'^0\^ < E(a;,w) . 
te[0,i] 

Combining this with (5.21) and the definition of I3, we finally conclude that (5.18) 
also holds for the term I3 (this time with a ~ 0), thus concluding the proof. □ 



Appendix A Some technical lemmas 

We give the proofs of some technical lemmas in this appendix. First, we show that 
Lemma 2.5 holds, which we reformulate here: 

Lemma A.l Let uj G W. Then, on any time interval bounded away from 0, the map 
t I— > Quit) is C°°. Furthermore, if we set f^it) = t-^Qojit), then we have /(^(O) — 
and for every T > 0, there exists a constant Mt such that II0.T.7 < -^^t||'^||(7.i5)- 

Proof. The proof follows from the calculations done in [HO07]. Since the function 
g(-) defined in (2.7) is smooth everywhere except at 0, it follows from (2.6) that mt is 
smooth on any interval bounded away from 0. Now a simple change of variables in 
(2.6) yields 

t'-^=,Hfltg'it).{-i)dr. (A.1) 

Since g'it) is smooth everywhere except at 0, for the map t 1-^ t ^it) to have finite 
7-H61der norm on [0, T], by (A.l) we need to only check the regularity of t g'it) at 
and at infinity. Indeed, if we show that the function t g'it) behaves like Oit) for i <C 1, 
and 0(t^^^/^) for < 3> 1, the conclusion follows from Lemma A.l and Proposition 
A.2 of [HO07]. We have, 

for a constant M. Notice that by Lemma A.l of [HO07] we have that git) behaves like 
Oit) for t < 1 and Oit"~^^'^) for t > 1. It can be seen easily that for i > 1, the term 

Jo ["+*')hIi/2 du behaves like Oit"'^/^). For i < 1, we have 
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< t u"-^/^ (l - (1 - uf'^-"^ du + M(t) (A.3) 

where M{t) is 0{t) for t < 1. Now since (1 - (1 - u)^/"^-") is 0{u) for u w 0, 
and diverges like 0((\ — u)^/^~^) for u « 1, the integral in (A.3) is finite. Thus from 
(A.2) it follows that tg'{t) is Oit) for t < 1 and we are done. □ 

Our next technical results yields a bound on the Holder constant of a function de- 
fined on a large interval as a function of the Holder constants of its restriction to smaller 
subintervals: 



Lemma A.2 For any N > I and any sequence of times Q ~ < ui < U2 < ■ ■ ■ < 

UN -I < um ~ T, the bound 

\\x\Wt^, < N'--" max (A.4) 

0<i<N — l 

holds. 



Proof. We proceed by induction. For = 1, the bound (A.4) holds trivially. Now 
assume (A.4) holds for some N G N and consider the sequence = ito < ui < U2 < 
• • • < utv-i < un < un+1 ~ T. If the supremum 5* = supj,_^(gjQ ^'^{fj'si^^^^^ is 
attained either for < s < ^ < or for wi < s < t < uat+i then we can use the 
induction hypothesis to conclude that 



||a;||o,T,7 = ll-j;||o,u„,7 V WAui 



< N'^ ^ max ||a;||„,,„ T.<(7V+1)^ ^ max ||a;||„,,„,, , 

and we are done. It remains to consider the case where the supremum S is attained for 
some s £ [0, ui] and t e [un, T]. We then have 

^ \x(ui) - x(s)\ + \x(U2) - x(ui)\ H h \x(t) - x(un)\ 

|t-s|7 

\ui - s\'^\\x\\q^uu'( \u2 - Ui\'''\\x\\m,u2n \t - Un\''\\x\\uj^,un^:^,j 

It-sl-y \t-s\f It-sl-y 

Since the function x H> is concave for 7 < 1 and since \ui — s\ + . . . + \t — u^l = 
\t — s\, we have 



< 



1 /|U1-S|'^ |m2-UiP , It-UATp 



< 



7V + 1 V \t-s\i \t-s\^ \t-s\-y 

1 |mi - s| 1 |m2-Ui| , 1 l^-UAflV 



.N + 1 \t-s\ N + 1 \t-s\ N+1 \t-s\ J (iV+1)'' 

Combining this with (A. 5), we obtain 

S<{N+1)^-^ max \\x\\^,.u,^,,^ , 

0<i<N + ' 

thus concluding the proof. □ 
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Finally, we provide the following interpolation result: 
Lemma A.3 The bound 



lo.T.oo < 2 max 



\T^^ '^II/IIl2[o,t], II/IIl2[o,t] "-/"IIo^t.^ 



holds for every 7 S (0,1] and every ^-Holder continuous function f : [0, T] i— >■ 

Proof Let /mm = infte[o,T] \f{t)\, so t) 
the definition of Holder continuity that 



Proof. Let ~ inffg[o,T] \f(t)\, so that < MJk^LZl, Therefore it follows from 



ll/lloXoo < /min + \\f\WT,^T'< < + ||/||o,T,^r^ . (A.6) 

For a > define the set 

L-f^{xe[0,T];\f(x)\>a}. 

11/11% 

By Chebyschev's inequality we have, for any a > 0, A(Lp < where A 

denotes the Lebesgue measure. Thus for any a > /min and x £ LJ, we have 

\f(x)\<a+\\f\\o,Tn\HL})\'' , 

so that 

||/||o,r,oo < a + a-'-'\\f\\o,T,4f\\l\o,T] ■ (A-7) 
Now we combine the estimates obtained in (A.6) and (A. 7). Making the choice a = 
\\f\&Jf\\l^,T] in (A-7) yields the bound 

||/||o,T,oo<2||/||o^^||/||jf|;^j, (A.8) 

provided that we have \\f\\^Jf\\]^^T] > > ./min- This is the case when 

ll./!lo,T,7 > \\f\\LHo.T]T-^-^ ■ When \\f\\o.T,^ < ||/|U2[o,T]r-^-i from (A.6) we 
obtain 

||/l!o,T,oo <2||/|U2[o,T]^-l/^ (A.9) 

thus completing the proof. □ 
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